Quantum path interference in the wavelength dependence of high-harmonic generation 
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We investigate the dependence of the intensity of radiation due to high-harmonic generation (HHG) 
as a function of the wavelength A of the fundamental driver field. Superimposed on a smooth power- 
law dependence observed previously we find surprisingly strong and rapid fluctuations on a fine A 
scale. We identify the origin of these fluctuations in terms of quantum path interferences with up 
to five returning orbits significantly contributing. 
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High harmonic generation (HHG) represents a versatile 
and highly successful avenue towards an ultrashort co- 
herent light source covering a wavelength range from the 
vacuum ultraviolet to the soft X-ray region [if. HHG 
has successfully opened new research areas, such as at- 
tosecond science 0, 0] and nonlinear optics in the XUV 
region 0, The fundamental wavelength A used in 
most of existing HHG experiments is in the near-visible 
range (- 800 nm). The cutoff law = \Ei\ + 3.17C/p, 
where \Ei \ denotes the binding energy of the target atom 
and Up = Fq/Auj"^ the ponderomotive energy {Fq: laser 
electric field strength), suggests that a longer fundamen- 
tal wavelength is advantageous to extend the cutoff to 
a higher photon energy, since Up quadratically increases 
with A. There is an increasing interest in the develop- 
ment of high-power mid-infrared (~ 2/Ltm) laser systems, 
e.g., based on optical parametric chirped pulse amplifica- 
tion. Along those lines the dependence of the HHG yield 
on A has become an issue of major interest. It has been 
commonly accepted that the spreading of the returning 
wavepacket would result in a A~'^ dependence of the HHG 
efficiency H as long as ground state depletion can be 
neglected fT]. Experimental findings Q have provided 
partial support. Recently, however, Tate et al. Q have 
reported a different wavelength-scaling of HHG between 
800 nm and 2 /im calculated with the time-dependent 
Schrodinger equation (TDSE) for Ar and a strong-field 
approximation (SFA) for He. They found a more rapidly 
decreasing HHG yield cx A~^' with 5 < x < 6. This sur- 
prising finding based on a somewhat limited number of 
data points motivated us to explore the A dependence in 
more detail employing two completely independent inte- 
gration methods of the TDSE to check for consistency 
and convergence. We have investigated the HHG for H 
and Ar on the level of single-atom response. Surprisingly, 
the harmonic yield does not smoothly decrease with fun- 
damental wavelength but exhibits rapid oscillation with 



a period of 6 — 20 nm depending on the wavelength re- 
gion. A semiclassical analysis based on the SFA reveals 
that the rapid oscillations are due to the interference of 
up to five different rescattering trajectories. Remarkably, 
averaged over the fast oscillations the smoothed yield fol- 
lows an approximate A~^ scaling, qualitatively consistent 
with the results of Tate et al. p]. 

We solve the atomic time-dependent Schrodinger equa- 
tion in a linearly polarized laser field in the length gauge. 
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where F{t) = Fq f (t) sm{Lut) denotes the laser electric 
field, f{t) is the envelope function and Vcs{r) the atomic 
potential. For H, Vcs{r) is the bare Coulomb potential 
while for Ar we employ a model potential within the 
single-active electron approximation which reproduces 
the binding energy to an accuracy of typically « 10^"^. 
We employ two complementary methods to solve Eq. ^ 
in order to establish reliable and consistent results. 

In the first method, Eq. ([T]) is numerically integrated 
using the alternating direction implicit (Peaceman- 
Rachford) method [llj with a uniform grid spacing Ar of 
6.25 X 10"^ a.u. In order to reduce the difference between 
the discretized and analytical wave function, we scale the 
Coulomb potential by a few percent at the first grid point 
The time step At is I/I6000 of an optical cycle for 
800 nm wavelength, i.e., 6.895 x IO~'^a.u.. This algorithm 
is accurate to the order of 0{At^). In the second method, 
the TDSE is integrated on a finite grid by means of the 
pseudo-spectral method [l3l| which is also accurate to the 
order of 0{At^). It allows for timesteps of the order of 
0.1 atomic units. The r coordinate is discretized within 
the interval [0, rmax] with a non-uniform mesh point dis- 
tribution. The innermost grid point is typically as small 
as 2.5 X 10"^ a.u., enabling an accurate description near 
the nucleus. A smooth cut-off function is multiplied at 
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each time-step to avoid spurious reflections at the border 
''max, while equivalently another cut-off function prevents 
reflections at the largest resolved energy i?max- Deeply 
bound, occupied states supported by the model potential 
are dynamically blocked during the time evolution 
We calculate the dipole acceleration d{t) = —df{z{t)), 
employing the Ehrenfest theorem through the relation 
d{t) = (?/'(r,t) I cos6l/r2 - F{t) \ ^{r,t))^, in which 
the second term can be dropped as it does not contribute 
to the HHG spectrum. 

For a direct comparison we adopt the laser parameters 
of with Ref. 0, with a flxed peak intensity of 1.6 x lO^"* 
W/cm^, a variation of A between 800 nm and 2 fim, and 
an envelope function f{t) corresponding to a 8-cycle flat- 
top sine pulse with a half-cycle turn-on and turn-off. We 
have checked that the fluctuations in the harmonic yield 
to be discussed below are not an artefact of this particular 
choice of f{t). They can be observed also for "smoother" 
pulse shape such as a sin^ pulse, provided that the pulse 
length is large enough to enable multiple returning tra- 
jectories. 

While we present the results obtained from the direct 
Fourier transform a{uj) of the dipole acceleration, we have 
confirmed the application of a Welch or Bartlett window 
[lit in the transformation hardly affects the results ex- 
cept for a constant factor. The HHG yield (defined as 
radiated energy per unit time, [lB|) integrated from 20 
to 50 eV, 
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calculated on a coarse mesh in A with a spacing of 
50 nm (Fig. [T]) falls off with a power law, AI oc A~^ 
(x « 4.8 — 5.5) for H and Ar, in qualitative agreement 
with Ref. . The two alternative integration algorithms 
employed in this work agree well with each other. Small 
discrepancies near 2 /xm are due to the difference in grid 
spacing and can be controlled by changes in the spacing 
near the origin. A power law (x « 5) results from the 
combination of two effects: the spreading of the return- 
ing wavepacket would give a: = 3 [6| for the overall yield. 
The increase of the cutoff Ec oc A^ results for a fixed 
energy interval (see Eq. ([2])) in an additional factor A~^. 

A closer look at Fig. [T] reveals the remarkable feature 
that the harmonic yield does not vary smoothly with A 
as may have been anticipated in the previous work, but 
strongly fluctuates. Slight change in fundamental wave- 
length may lead to variations of the yield by a factor of 
2 to 6. Such rapid fluctuations imply that a reliable A 
dependence can only be established by employing a much 
flner A grid. Moreover, the notion of a simple power law 
scaling itself is called into question and can apply, if at 
all, only after averaging over fluctuations. The fluctua- 
tions are not specific to hydrogen but appear for argon 
(Fig. IHb)) as weU. 
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FIG. 1: Integrated harmonic yield A J between 20 and 50 eV 
as a function of A calculated on a coarse mesh with AA = 
50 nm. •: Peaceman-Rachford method; □: pseudo-spectral 
method, solid line: fit AI oc A~^: a) hydrogen, b) argon. 
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FIG. 2: (Color online) Fluctuations of the harmonic yield AI 
as a function of the fundamental wavelength A and the photon 
energy E for hydrogen. An oscillatory behaviour of the yield 
as a function of A emerges after integrating over vertical strips 
(upper panel and Fig. |3ja)). 



FigurelUpresents a two-dimensional zoom into the fine- 
scale variations calculated on a mesh with AA = 1 nm, 
as a function of A and the photon energy E. Fluctua- 
tions appear which form vertical "ridges" spanning sev- 
eral harmonic orders. Consequently, after integrating 
over vertical strips we arrive at a fiuctuating function 
A/(A) (upper panel). Such variations in quantities av- 
eraged over a large energy interval are expected when 
few (quantum) orbits with a well-defined time structure 
significantly influence the spectrum ( "periodic orbit spec- 
tra" ) [l3| • A superposition of oscillations of different fre- 
quencies emerges (Fig. [3]), the dominant of which has a 
fluctuation scale SX of « 20 nm near a driver wavelength 
of A = 1000 nm, and (5A ~ 6 nm near A = 2000 nm. 
Similar oscillation patterns can be observed for argon on 
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Fig. [Sljc). Apparently, they are largely independent of 
the atomic species. Such oscillations are obviously the 
result of interference effects. 

Many features of HHG can be intuitively and even 
quantitatively explained in terms of quantum trajecto- 
ries 18 I which represent the semiclassical three step 
model 19|. The main contribution to the HHG spec- 
trum comes from those electronic quantum paths that 
correspond to classical returning trajectories ionized at a 
certain time ti and recombining with the parent ion at a 
later time tf. In order to identify the origin of the inter- 
ference structures we apply a semiclassical model based 
on the strong field approximation (SFA) . In this 

model, the time-dependent dipole moment d{t) can be 
expressed as [ISl 
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i.e. a sum over paths P that start at the moment of tunnel 
ionization ti with amplitude 6ion(ii)j evolve in the laser 
field - e~*'^^^*''*^^ - and recombine upon rescattering at 
the core at time tf with the amplitude Crec(i/)- The sum 
over all possible electron trajectories recoUiding at time 
tf may be large, but in practice is limited by wave-packet 
spreading. We consider up to 16 possible ionization times 
ti for each individual tf. 

Interference oscillations are controlled by the evolu- 
tion phase, the semiclassical action of the path P, which 
reads: 



Sp{ti,tf) 



ip + A{t')f 



■dt' + g-\E,\{tf -t,) (4) 



\Ei\ is the ionization potential (binding energy) of the 
atom and A{t) the laser vector potential defined by 
A{t) = — F(t')dt' . p is the classical momentum of the 
returning trajectory, g = 1.3 is a correction factor that 
accounts for the dependence on the tunneling time, modi- 
fying the energy of the recoUiding photons [a, [21 1 . When 
including up to five returning paths, the semiclassical 
calculation can reproduce the modulation depth, mod- 
ulation frequency, and the approximate phase of the A 
oscillations reasonably well, thus unambiguously estab- 
lishing the quantum path interference as the origin of 
the fluctuations (Fig. [S^b)). Consideration of additional 
trajectories leads only to minor modifications with no 
qualitative differences. Setting g = 1.0 yields only a 
phase shift of the predicted oscillations while retaining 
the modulation frequency and depth. 

Very recently, first experimental evidence for the inter- 
ference between the short and long paths was presented 
[2^ in the intensity dependence of the HHG yield. Re- 
markably, for the present A dependence the frequently 
discussed short and long trajectories (dotted curve in Fig. 
[3]Jb)) are insufficient to account for the interference os- 
cillations. Convergence is approximately reached only 
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FIG. 3: Variations in the integrated harmonic yield (20 to 50 
eV) in a narrow range of A. a) comparison between the TDSE 
solution with the SFA including 5 quantum paths for hydro- 
gen, b) build-up of the interference pattern with increasing 
number of quantum trajectories within the SFA. Thin dashed; 
1 (short) trajectory. Thin dotted: 2 (short and long) trajec- 
tories, Bold dashed: inclusion of the shortest 5 trajectories, 
c) as a) but for argon. 



when at least five trajectories are included. Moreover, 
the presence of the oscillations is independent of the par- 
ticular choice of the envelope function f{t). Note that the 
Fourier broadening of the few-cycle driving field exceeds, 
on a wavelength scale, the period SX of the modulation. 
This somewhat surprising finding is a direct consequence 
of the quantum path interference. As long as the few- 
cycle pulse permits the generation of a set of a few quan- 
tum paths in subsequent half-cycles, the overall temporal 
characteristics of the driver pulse is of minor importance, 
though the latter will infiuence the detailed shape of the 
interference pattern. Preliminary calculations for pulse 
propagation in one dimension, accounting for the geomet- 
ric Guoy phase, show that inteference oscillations persist 
in loose focus geometry. 

The principal modulation length of the harmonic yield 
SX is a function of the wavelength A itself. It approaches 
w 6 nm near a wavelength of 2 fim (Fig. 3]). For a sim- 
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is different from the generally accepted A ^ scaling, but 
is close to that reported in Ref. [§]. 
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FIG. 4: Variation of the modulation period SX with the fun- 
damental wavelength A for hydrogen. ■: TDSE, SPA. 
Solid line: A~^ scaling (see Eq. (|6])). 



pie estimate for the scaling of SX with A we note that the 
semiclassical action in Eq. ^ has its largest contribution 
from the A{t)'^ term in the strong field case. Hence, the 
phase difference between the shortest and longer tajecto- 
ries due to the semiclassical action can be approximated 
by Sp K, Up ■ Tf where 'Tj_— tf — is the flight time of 
the electron trajectory [l8|]. Bp scales approximately as 
A'^. The period of the modulation corresponds to a phase 
change of Sp by 2tt. Accordingly, 
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This estimate can be improved when using the full ex- 
pression for Sp; nevertheless even Eq. © predicts the 
A dependence of the modulation length remarkably well 
(see Fig. SI). 

In conclusion, we have found that the fundamental 
wavelength dependence of HHG in the single-atom re- 
sponse features surprisingly strong oscillations on fine 
wavelength scales with modulation periods as small as 
6 nm in the mid-infrared regime near A = 2 /im. These 
oscillations have been established for both hydrogen and 
rare gas targets (as an example argon is shown in this 
work) using two complementary integration algorithms 
for the TDSE. Thus, even a slight change in fundamen- 
tal wavelength leads to variations in the HHG yield by a 
factor of up to ~ 6. According to our semiclassical anal- 
ysis based on the SEA, this unexpectedly rapid variation 
on a fine scale is the consequence of the interference of 
different rescattering trajectories. We have to take ac- 
count of up to five returns to reproduce the results of the 
quantum simulations. This confirms the significance of 
higher-order returns of the electron wavepacket [9| . On a 
large A scale, apart from the rapid oscillation, our TDSE 
results show that the HHG yield at constant intensity de- 
creases as A~^ with a; w 5 for H and Ar. This dependence 
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